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ABSTRACT 



Aims. We investigate the effects of coherent isotropic continuum scattering on the formation of spectral lines in local thermodynamic 
equilibrium (UTE) using 3D hydrodynamical and ID hydrostatic model atmospheres of red giant stars. 

Methods. Detailed radiative transfer with coherent and isotropic continuum scattering is computed for 3D hydrodynamical and ID 
hydrostatic models of late-type stellar atmospheres using the SCATE code. Opacities are computed in UTE, while a coherent and 
isotropic scattering term is added to the continuum source function. We investigate the effects of scattering by comparing continuum 
flux levels, spectral line profiles and curves of growth for different species with calculations that treat scattering as absorption. 
Results. Rayleigh scattering is the dominant source of scattering opacity in the continuum of red giant stars. Photons may escape from 
deeper, hotter layers through scattering, resulting in significantly higher continuum flux levels beneath a wavelength of A < 5000 A. 
f~| , The magnitude of the effect is determined by the importance of scattering opacity with respect to absorption opacity; we observe the 

^ |i largest changes in continuum flux at the shortest wavelengths and lowest metallicities; intergranular lanes of 3D models are more 

I ' strongly affected than granules. Continuum scattering acts to increase the profile depth of LTE lines: continua gain more brightness 

O ' than line cores due to their larger thermalization depth in hotter layers. We thus observe the strongest changes in line depth for high- 

excitation species and ionized species, which contribute significantly to photon thermalization through their absorption opacity near 
C/3 , the continuum optical surface. Scattering desaturates the line profiles, leading to larger abundance corrections for stronger lines, which 

^ ' reach -0.5 dex at 3000 A for Fell lines in 3D with excitation potential = 2eV at [Fe/H] = -3.0. The corrections are less severe 

for low-excitation lines, longer wavelengths, and higher metallicity. Velocity fields increase the effects of scattering by separating 
emission from granules and intergranular lanes in wavelength. ID calculations exhibit similar scattering abundance corrections for 
^ . weak lines, but those for strong lines are generally smaller compared to 3D models and depend on the choice of microturbulence. 

' Conclusions. Continuum scattering should be taken into account for computing realistic spectral line profiles at wavelengths A < 

^ 4000 A in metal-poor giant stars. Profile shapes are strongly affected by velocity fields and horizontal inhomogeneities, requiring a 

. treatment based on 3D hydrodynamical rather than classical ID hydrostatic model atmospheres. 
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1. Introduction local coupling enables the approximation of local thermody- 

,,. _ .. . namic equilibrium (LTE), which allows a tremendous simplifi- 

Spectral hne formation is an important discipline in the field of ^^^.^^ ^.^^ formation computations and is therefore frequendy 

stellai- astronomy. It is used as a diagnostic tool for numerical ^^^^ However, LTE is in general a bad approximation, and it can 

^ models of stellar atmospheres by comparing their predictions ^^^j j^^^ .^^ ^^^^^ ^^.^^ .^^ ^^j.^j j^^^ ^^^^ 

H , with observations, as well as for measuring quantities of as- .^^^^^ ^j^^^^ ^^^j ^^^^^^^^^ quantities ai-e needed. A more 

. ^ . trophysical interest such as stellar parameters, chemical abun- ^^^^^^^^ j^^^^.^te treatment involves detailed atomic models, 

dances, surface velocity fields and many more Specti-al lines ^j^^^j^ ^ ^ significant amount of additional physical data 

sample the physical conditions m stellar atmospheres in a wide ^^^^^ discussion in&nd 2005). Owing to its much 

height range and can thereby reveal useful information about the j^^^^. complexity and computational demands, 3D LTE line for- 

atmospheric structure. mation is nevertheless a useful tool for cases that are known to 

The physics of line formation involves many diff^erent pro- ^j^^^ ^^^^ vanishing departures from LTE, such as the often- 

cesses that require a detailed treatment of the interaction be- , „ , . , , t^tt i- ^Tr,r, * 

J- J „ r. J .u . ir .u used forbidden [0 1] hne at 6300 A. 
tween radiation and matter. Beside the transition itself, these 

processes include, e.g., Doppler-shifts through thermal and non- Late-type stellar atmospheres are bounded by the underly- 

thermal particle motion, collisional interaction, photoionization ing convection zone, which is responsible for horizontal inho- 

etc. The physical state of matter depends in general not only mogeneities in gas temperature and pressure, and flow fields 

on local conditions, such as the gas temperature and pressure, that affect the line p rofile shape through Doppler-shifts (e.g., 

but also on conditions in other parts of the atmosphere that INordlund et alj 120091) . Accurate line formation computations 

are connected through the radiation field. Ignoring such non- need to take this 3D structure of the atmosphere into account to 
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obtain correct line profiles. The complexity of 3D calculations is 
rewarded by the elimination of additional parameters such as mi- 
croturbulence and macroturbulence (e.g. Asplund et al. 2000b). 

Metal-poor halo stars are particularly interesting targets 
for stellar abundance analyses: they are among the oldest ob- 
jects in the Galaxy and largely preserve the chemical compo- 
sition of the ISM in their envelopes, allowing detailed stud- 
ies of nucle osynthesis processe s and Galactic chemical evo- 
lution (e.g., iBeers & Christhel?! 2005). Most abundance stud- 
ies of metal-poor stars rely either on ID LTE or on ID non- 
LTE syntheses. The effects of 3D lin e formation h a ve no t 
yet been fti lly explore d, although, e .g., Asplund et al] (1 19991) . 
ICollet et al.1 (12007) and lBehara et aH (12010) ' have shown that the 
temperature stratification and inhomogeneous structure of 3D 
model atmospheres have an important effect on line formation 
in metal-poor dwarf stars and giant stars. Molecular species, 
low-excitation lines and atomic minority species are thus par- 
ticularly affected by 3D effects through their high temperature- 
sensitivity; molecules are important instruments for determin- 
ing CNO abundances of such stars. Full 3D non-LTE studies of 
lithium abundances i i i metal-poor stars have been presen ted by 
lAspl und et al.V2003VBarkl em et al.l(l2003l) .ICavrel et al.'(2007) 
andjSbordone et al., (.2010). fo llowing earlier work of .Kiselmam 
(Il997h and lUitenbroekl (119981) for the Sun. 

An important issue for the accuracy of abundance analyses of 
metal-poor giants is the role of Rayleigh scattering opacity for 
continuum formation in the blue and UV wavelength regions, 
where transitions of many important elements are found. Such 
scattering processes can increase the thermalization depth of 
photons, which may escape from deeper, hotter layers in the star, 
resulting in higher continuum flux. Treating continuum scatter- 
ing as absorption instead of computing radiative transfer with 
scattering in ID chemical abundance analyses leads to signifi- 
cantly larger abundances and trends with transition wavelength 
(ICavrel et al.l l2004'. Bi hain etaDl2004t iLaTetaT 2008). 

We investigate the effects of coherent scattering in the back- 
ground continuum by computing 3D LTE line profiles for a se- 
lection of typical transitions, including neutral and singly ion- 
ized iron, as well as various molecular species. Section |2] dis- 
cusses important aspects of the underlying radiative transfer 
model for spectral line formation, which builds on 3D time- 
dependent radiation-hydrodynamical model atmospheres of red 
giant stars discussed in Sect.|3] We analyze continuum formation 
with scattering in Sect. |5] investigate its effects on 3D LTE line 
formation in Sect. |6]and derive 3D-3D abundance corrections 
for typical spectral lines in Sect. |7] 3D-1D abundance correc- 
tions computed with continuum scattering and treating scattering 
as absorption are compared in Sect.|8j we conclude in Sect.|9] 

2. Line formation with continuum scattering 

We solve the 3D time-independent radiative transfer problem for 
stellar spectral lines and a background continuum with coherent 
isotropic scattering. Excitation level populations, ionization and 
molecule formation needed to derive gas opacities are computed 
assuming local thermodynamic equilibrium (LTE); the line pro- 
files include collisional broadening by neutral hydrogen and 
Doppler-shifts through macroscopic velocity fields (see Sect. |4] 
and Appendix lAlfor a description of the numerical methods). 

The time-independent radiative transfer equation evaluates 
the monochromatic specific intensity ly in direction h at fre- 
quency V in the observer's frame: 

— = -/v(Tv) + Sy(Ty), (1) 

dTy 



where 5v is the monochromatic source function. The optical 
depth dTy of a photon path length ds on a light ray is defined 



dTy = {xl+x')ds, (2) 

with the line opacity Xv the continuum opacity x'^- The 
source function Sy at optical depth Ty includes contributions 
from spectral lines and continuum processes: 



Sy^ 



Xy 



-B + 



X'y+X' 



X'y+X' 



■([l-e=]7, + e^B); 



(3) 



B is the Planck function, e'^ is the continuum photon destruction 
probability and Jy is the monochromatic mean intensity. Note 
that frequency subscripts have been dropped in Eq. ([T]i through 
Eq. (O for all quantities that depend only weakly on frequency 
across a line profile and are therefore assumed constant in our 
computation. 

The line opacity x\, may include contributions from several 
transitions of different species. At each point in the atmosphere, 
the individual line profiles are evaluated in the local rest frame 
of the gas around the center frequency of the transition. In 
the observer's frame, the profiles appear Doppler-shifted through 
the non-relativistic velocity field u predicted by the 3D hydrody- 
namical model atmosphere. The line center frequency vq is then 
given by 



(4) 



vo = 1 + 



in the observer's frame f or a light ray in direction h (see, e.g., 
Mi halas & Mihalas ''1984); c is the speed of light. Relative to the 
frequency v of photons that stream in the same direction as the 
flowing gas, the absorption profile appears blue-shifted; the same 
holds for line emission from the observer's point of view. In the 
solar atmosphere, upflowing gas in the bright granules dominates 
the emitted radiative flux (S tein & No rdlund 1998), and the re- 
sulting Doppler-shifts along the line of sight distort profile bisec- 
tors dominantly to wards higher frequenc ies (see the discussion 
in Sect.|6]and, e.g.. 'As plund et al.ll2000bl) . 

The opacity quotients x\l{x\+x'') and x''l{x\+x'') in 
Eq. ^ yield the probabilities that a photon was created by 
a line transition or by a continuous process. The direction- 
dependence of line opacity in the observer's frame through 
Doppler-shifts (Eq. (|4]i) induces a direction-dependence of these 
emission probabilities, and the combined source function S y be- 
comes anisotropic. 

The continuum photon destruction probability e'^ is defined 
by the ratio 



e = 



K- + cr- 



(5) 



with the continuous absorption coefficient k'^ and the continu- 
ous scattering coefficient cr'^. In analogy to the above described 
opacity quotients, e'^ can be interpreted as the probability that 
a photon was fed into the beam from the thermal pool rather 
than from a scattering event. In late-type stellar atmospheres, 
Rayleigh scattering on H I atoms and electron scattering mainly 
contribute to continuous scattering opacity, while continuous 
absorption includes many bound-free and free-free transitions 
of hydrogen, helium and metals (see, e.g., the discussion in 
Hayek et al. 2010). 

Continuous scattering may be treated as a coherent mech- 
anism to very good approximation when the radiation field Jy 
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varies only slowly with frequency. In the presence of spec- 
tral lines, Jy exhibits much stronger frequency-dependence, 
and scattering leads to a complex coupling of the radiation 
field in freque ncy and angle t hrough Dopple r - shifts (see, e.g., 
discussions in iMihalasI Wf^ iPeraiahl l200lh . lAuer & MihalasI 
(Il96 8a.b) analyzed the effects of non-coherent electron scatter- 
ing on line formation in early-type stars. They assumed redis- 
tribution through thermal motion of the electrons, which occurs 
across a very wide frequency band compared to thermal redistri- 
bution in the line: the ratio between the thermal profile widths of 
electrons with mass and line-forming atoms with mass matom 
is given by Vwiatom/'We ~ 42.7 VA, where A is the atomic mass 
number. In late-type stars, Rayleigh scattering dominates scat- 
tering opacity at the continuum optical surface in the blue and 
UV wavelength regions, and the ratio of Doppler widths scales 
only with a/A < 15. Although electron scatterin g is an important 
continu ous opacity source in early-type stars, lAuer & MihalasI 
(Il968hh find a relatively small effect of non-coherence on the 
wings of Hell profiles. We therefore expect the assumption of 
coherent scattering to be a reasonable approximation for the 
rather weakly scattering background continua of red giant stars. 

Rayleigh scattering and electron scattering cross-sections 
depend on the scattering angle through the angular redistri- 
bution function (In- cos ^ 6), which ma y be neglected to good 
approximation (see, e.g.. iMUiala J 1 978h . 



3. 3D radiation-hydrodynamical model atmospheres 

Various studies have shown that the formation of spectral lines in 
the atmospheres of late-type stars can be severely affected by the 
presence of inhomoger ieities in t h e teni p erature structure and 
veloc ity fields (e.g., 'Nordlund' '1980^, 'Dravins & Nordlund' 
199Ct [A splund et al. 2000b; Steffen & Holweger 2002; 
AUende Prieto et al.1 120021: ICaffau et al.l l2008l) ." Synthetic 
line formation computations based on the current generation of 
3D radiation-hydrodynamical models are capable of providing 
realistic predictions of the observations; see Perei ra et al.l 
(2009a.b) for the important case of the Sun. Classical ID 
hydrostatic model atmospheres, which are still widely used, 
simply cannot achieve this degree of realism. 

An important aspect in this context is the coupling of 
radiative transfer and hydrodynamics in 3D models, which 
leads to significantly cooler gas temperatures above the sur- 
face ofmetal-goor stare than predicted by ID hydrostatic mod- 
els (lAsplund et al.|[T999l) . Low metallicity strongly reduces ra- 
diative heating in these objects, and the photospheric tem- 
perature stratification steepens towards an adiabatic gradient, 
while the assumption of radiative equilibrium that is inherent 
to ID models keeps the gradient shallow. The lower tempera- 
tures of 3D models have a strong impact on the formation of 
molecules and thu s ori th e predicted molec ular line strengths 
lAsplund & Garcia PerezllSoOll: l(5olM et^^ allllOOTl: iBehara et al l 
120101) : molecules are frequently used for de termining th e abun- 
dances of carbon, nitrogen and oxygen (e.g.. lSpite et al.l l2006). 

We base our analysis on 3D radiation-hydrodynamical model 
atmospheres of red giant stars with similar effective tempera- 
tures, a surface gravity of log^ = 2.2 (in cgs units) and metal- 
Ucitie^] between [Fe/H] = -3.0 and [Fe/H] = 0.0 (see Tabled; 
the chemical abundance mixture adopts the solar composition of 



' Metallicity is defined here in spectroscopic notation as the loga- 
rithmic iron abundance relative to the Sun, [Fe/H] = logmCA'pe/A^n). - 

logioCiVpe/A'H)©. 



lAsplund et al.l (l2009l) . scaled by metallicity and with a -1-0.4 dex 
enhancement of a-elements for the most metal-poor model. 

The simulations were c reated with the StaggerCode 
(iNordlund & Galsgaardll 19951) . which solves the coupled equa- 
tions of compressible hydrodynamics and time-independent 
radiative transfer, computing radiative heating rates in LTE. 
Line-blanketing is approximated using the opacity binning 
method (Nordlund 1982; SkartHen 2000) with 4 opacity bins. 
Collet et al. ( 201 li) demonstrate that the effects of coherent con- 
tinuum scattering on the temperature structure of metal-poor gi- 
ant stars with [Fe/H] = -3.0 and [Fe/H] = -2.0 may be repro- 
duced to good approximation by removing the scattering contri- 
bution from the opacities in the optically thin parts of the atmo- 
sphere and using a Planck source function, which we applied to 
the models used for the analysis to reduce computational load. 
The deviations in average temperature from simulations with co- 
herent scattering radiative transfer remain below < 50 K in the 
photospheres of the two most metal-poor giants; see lCollet et al.l 
(1201 Ih for further details. 

The effective temperature of 3D radiation-hydrodynamical 
model atmospheres is an observable rather than a parameter: it 
is determined by the entropy of the inflowing gas at the bottom 
of the simulation, but exhibits some temporal variation due to the 
convective motions. The time-averaged (Teff) therefore varies 
between the different model atmospheres, which is not critical 
for our differential investigation of continuum scattering effects. 
After the initial scaling of the model atmosphere to reach the 
desired stellar parameters, each simulation was run until the at- 
mospheric stratification reached a quasi-steady state with little 
temporal variation. 

Figure [T] shows the temperature distribution of an arbitrary 
snapshot of the [Fe/H] = -3.0 model as a function of optical 
depth at 5000 A. The atmospheres cover several pressure scale 
heights above and below the continuum optical surface of the 
stars to include the relevant line-forming regions and to avoid 
boundary effects on the granulation flow. The simulation do- 
mains have a resolution of 240x240x230 grid points and assume 
periodic horizontal boundaries. Synthetic curves of growth were 
computed using time-series of simulation snapshots that span 
several periods of the fundamental pressure oscillation mode; a 
representative sample of line profiles is thereby obtained that in- 
cludes statistical variation in the predicted line strengths through 
convective motions. 

Rece nt calculations of OH line profiles by 
Gonzale z Hernandez et al.l (2010) based on 3D model at- 
mospheres and radiative transfer with 12 opacity bins yield 
smaller 3D-1D LTE abundance coiTections compared to their 4 
bin models. The refined opacity treatment leads to a shallower 
temperature gradient, which reduces molecular equilibrium 
populations in the atmosphere. The OH abundances produced 
by the 12 bin 3D models are thus closer to those derived 
from ID flux-equilibrium models. We test the importance 
of the opacity treatment for the most me tal-poor giant sta r 
where 3D abundance effects are largest (see ICollet et al.ll2007h 
using a 3D model constructed with a 12 bin opacity table. 
The opacities are sorted by Rosseland optical depth of the 
monochromatic opt ical su rfaces and by wavelength, similar to 
iGonzalez Hernandez et al.l (l2010i) . The time-averaged effective 
temperature of the model is {T^s) - 5058 K and thus very close 
to the 4 bin model (see Table[T]). 

We also include MARCS model atmospheres with the same 
stellar parameters as the time-averaged 3D models (see Table 
[TJ in our analysis to compare the effects of scattering between 
3D hydrodynamical and ID hydrostatic atmospheres. Each ID 
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Table 1. Stellar parameters of the 3D radiation-hydrodynamical model atmospheres. 



(T-eff) [K] 


log^ [cgs] 


[Fe/H] 


[am 


5100 


2.2 


-3.0 


+0.4 


5051 


2.2 


-2.0 


0.0 


4730 


2.2 


-1.0 


0.0 


5063 


2.2 


0.0 


0.0 



model is converted into a 3D box with zero gas velocity, hor- 
izontally homogeneous layers and with the same T-tsooo rela- 
tion; line formation is computed using the exact same method as 
in the 3D case. 




10"'* 10"^ 10° 10^ 

■''sooo 

Fig. 1. Temperature distribution as a function of vertical contin- 
uum optical depth r at 5000 A, computed for an arbitrary snap- 
shot of a 3D radiation-hydrodynamical model atmosphere with 
<reff) = 5100K, logg = 2.2 (cgs) and [Fe/H] = -3.0. The 
dashed line indicates the average temperature at each optical 
depth. 



4. 3D line formation computations 

Synthetic continuum flux and line flux profiles are calcu- 
lated using the SCATE code. The program first computes 
the monochromatic continuum source function with coherent 
isotropic scattering using a short characteristics-based radiative 
transfer solver with Gauss-Seidel-type approximate A-iteration 
dTruiillo Bueno & F abiani Bendicho 1995). A detailed descrip- 
tion of the implem entation can be found in Appendix [Aland in 
iHavek et al.l (l2010l ). Continuum opacities and photon destruc- 
tion probabilities are looked up in precomputed tables, line opac- 
ities are calculated during runtime. All quantities assume the 
LTE approximation. 

The numerical method takes Doppler-shifts of line profiles 
into account, which influence absorption and line-to-continuum 
photon emission probabilities as a function of ray direction, 
gas velocity and frequency according to Eq. Q. Including this 
inherent coupling is essential for correctly predicting the im- 
pact of scattering on the local radiation field and consequently 
on the profile shapes (see Sect. |6]l. A minimum resolution in 
solid angle and frequency is furthermore needed to reproduce 
the effects of Doppler-shifts. For the case of continuum scatter- 
ing with LTE lines in late-type stellar atmospheres, where gas 



flow reaches only moderate velocities, ICarlsonI (1 19631) quadra- 
ture with 24 ray directions provides sufficient accuracy; comput- 
ing radiative transfer with 48 angles changes the spatially av- 
eraged flux level by < 0.3%. Carlson quadrature also has the 
advantage of rotational invariance, which avoids directional bias 
for determining the local mean radiation field Jy. Line profiles 
for synthetic curves of growth are computed with typically 40 
frequency points. We use linear (non-logarithmic) interpolation 
of quantities between the grid of the model atmosphere and the 
characteristics grid, as well as for the source function integral in 
the optically thin regime. Second-order interpolation is applied 
to the source function integral in optically thick regions to cor- 
rectly recover the diffusion approximation. 

Once the continuum source function has converged, it is 
passed to a second radiative transfer solver, which computes 
outgoing specific intensities on characteristics that span across 
the entire 3D atmosphere cube. This method has the advantage 
of reduced numerical diffusion with respect to the short char- 
acteristics method when angle-resolved surface intensities are 
needed. Local cubic logarithmic interpolation translates the rel- 
evant quantities onto the tilted characteristics grid. The radiative 
transfer equati on is then solved along vertical columns using the 
lFeautriei](ll96 4) method to obtain surface intensities. 

The grid resolution of the model atmosphere is an impor- 
ta nt issue for 3D line fo rmation calculations (see the discussion 
in lAsplund et al.ll2000al) . The axis spacing in the vertical dom- 
inates the accuracy of the solution, due to the strong tempera- 
ture gradients near the stellar surface. Our calculations are based 
on a mesh with 120 x 120 x 230 grid points, and we refine the 
vertical grid by automatic insertion of additional layers where 
needed to obtain robust intensity and flux profiles. Computation 
of the scattering source function excludes the deepest layers of 
the 3D model atmosphere, which are optically thick and domi- 
nated by local thermal radiation. The Feautrier-type solver starts 
at a continuum optical depth logjQ t*^ » 2 on each individual ray 
to integrate the entire contribution function of continuum and 
line emission. 

Stellar spectroscopy measures the radiative flux integrated 
over the stellar disk in most cases of interest. The monochro- 
matic flux that leaves the surface on the visible hemisphere 
into the observer's direction is given by 

Fv,z= {Iy)iti,cf>)ndfid<p, (6) 

Jfi=0 J</>=0 

where jj = cos Q is the projection factor, Q is the polar angle, ^ is 
the azimuth angle and (/y) is the horizontal average surface in- 
tensity of the simulation cube. The 3D hydrodynamical model 
atmospheres are interpreted as local statistical representations 
of stellar surface convection. Each angle pair (yu, 0) for which 
radiative transfer is computed in the 3D cube corresponds to a 
different position on the disk. The flux integral (Eq. (|6]l) is thus 
equivalent to an integral over the stellar surface as seen by the 
observer; Fy.z thereby automatically includes the limb darkening 
effect. Including stellar rotation in 3D spectral line formation 
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computations requires additional consideration (lLudwigll2007h . 
we therefore assume zero rotation (v sin / = 0) for simplicity. 

We approximate the flux integral (Eq. (|6]l) using Gauss- 
Legendre quadrature for the polar angle and the trapezoid rule 
for the azimuth angle; 4 x4 ray angles reproduce the surface flux 
with good accuracy. Doubling the number of polar angles or az- 
imuth angles changes the spatial average flux level by < 0.4 %. 

5. The effects of scattering on the continuum flux 

We first seek to determine the wavelength range in which con- 
tinuum scattering contributes sufficient opacity to influence con- 
tinuum flux levels. Radiative transfer is computed in 3D with 
scattering opacity treated as absorption and as coherent scatter- 
ing for a set of wavelength points with logarithmic distribution 
in the range 1500 A < A < 8000 A, taking only continuum opac- 
ity into account; Fig.|2]shows the resulting flux distributions for 
the different model atmospheres (see Table [T] for stellar param- 
eters). It is clear that continuum scattering is only significant in 
the UV and blue regions in all cases, where Rayleigh scattering 
on HI atoms is mostly responsible for the increased flux lev- 
els compared with the case of treating scattering as absorption; 
the /I "^-dependence of the Rayleigh cross-section and thermaliz- 
ing bound-free transitions of primarily various metals at shorter 
wavelengths limit the spectral range. Scattering is therefore com- 
pletely negligible in the infrared. 

Continuum scattering evidently increases the thermalization 
depth in the UV, as a significant fraction (1 - e^) of outward 
streaming photons from deeper, hotter layers is scattered instead 
of absorbed. Figure |3] shows the vertical variation of photon de- 
struction probabilities (e'^) in arbitrary snapshots of the model 
atmospheres, averaged over surfaces of constant vertical contin- 
uum optical depth t\ at the three wavelength points indicated 
by vertical dotted lines in Fig.|2] For the [Fe/H] = -3.0 model 
(upper left panel of Fig. O, all photons that interact with the 
gas above optical depth rjj > 10 are absorbed and thermalize. 
At 3000 A, the atmosphere becomes slightly translucent before 
the optical surface (rjj = 1) is reached, allowing photons to 
escape from larger depths than in the case where scattering is 
treated as absorption, while photons at 4000 A and 5000 A are 
still trapped. Higher up in the atmosphere, scattering becomes 
increasingly important at all wavelengths. Photon destruction 
probabilities grow with metallicity, but exhibit only weak vari- 
ation at the optical surface for the models with [Fe/H] < -1.0. 
The slightly lower effective temperature of the [Fe/H] = -1.0 
model reduces ^ above the surface relative to the other model 
atmospheres due to the temperature-dependence of absorption 
opacity (see the discussion below). At [Fe/H] - 0.0, absorption 
dominates almost everywhere in the photosphere and beneath. 
The H I (« = 2) bound-free absorption edge at 3647 A (the so- 
called Balmer jump) causes a depression in the continuum flux 
on the blue side. Absorption through photoionization increases 
the photon destruction probability and moves the thermalization 
depth outward into cooler layers, causing a spike in the flux de- 
viation A (f^) (upper panels of the plots in Fig. |2]l. 

The granulation flow in 3D radiation-hydrodynamical simu- 
lations produces strong horizontal variations of the outward ra- 
diative intensities between the hot, bright granules and the cool, 
dark intergranular lanes. The left panel in Fig. |4] shows contin- 
uum intensities in the stellar disk center at 3000 A, computed 
with scattering as absorption for an arbitrary snapshot of the 
[Fe/H] = -3.0 model. The shapes of granules and intergranular 
lanes are similar to solar granulation, but extend to much larger 



spatial scales (ICoUet et al.l2007h . Note that radiative emission in 
the dark intergranular lanes is comparatively small but nonzero, 
which is important for the formation of line profiles (see the dis- 
cussion in Sect. |6]l. The right panel in Fig. |4] shows disk center 
continuum intensities computed with scattering; both panels in 
Fig. |4] are normalized to the same average intensity. While the 
overall morphology of the granulation pattern is almost iden- 
tical, the surface intensities appear slightly brighter due to the 
larger thermalization depth of the coherent scattering case. 

The left panel of Fig. |5] quantifies the spatially resolved in- 
tensity ratio of the scattering and absorption calculations, show- 
ing the correlation with continuum surface intensities; the darker 
intergranular lanes gain more brightness in proportion to their in- 
tensity than the brighter granules. The reason for this variation 
is the temperature-dependence of the continuum photon destruc- 
tion probabilities e'^, visible in their correlation with continuum 
intensity at local optical surfaces (t'^ - 1) in each column of 
the model atmosphere (right panel of Fig. |5j: the temperature- 
dependence of Rayleigh scattering opacity is weak if the gas 
is cool enough that hydrogen is dominantly neutral, cr is there- 
fore very similar at continuum optical surfaces in granules and 
lanes. Contrary to that, thermalizing opacity varies strongly with 
temperature. Absorption is thus dominant in hot granules, while 
scattering opacity is important in lanes. The result is a larger 
thermalization depth with respect to the optical surfaces in the 
lanes compared to granules and a relatively stronger intensity 
gain through scattering. The overall contrast in the continuum 
surface intensities reduces through scattering: while the horizon- 
tal average continuum intensity (F) at 3000 A of the snapshot 
shown in Fig. |4] increases by s; 27 % with respect to treating 
scattering as absorption, the relative rms variation at this wave- 
length decreases from I^^^J (F) = 0.51 to I^^J (F) = 0.43 when 
coherent scattering is included in the radiative transfer computa- 
tion. 



6. The effects of scattering on spectral lines in 3D 

We now add spectral lines from different atomic and molecu- 
lar species to the continuum opacities. It is sufficient to restrict 
the wavelength range to 3000 A < A < 5000 A, in which ra- 
diative flux is accessible for measurements with current instru- 
ments and in which continuum scattering is important, as it was 
demonstrated in Sect.|5] The profil e bro adening in all Fe line for- 
mation computations assumes the lUns51dl ( 11955 ) approximation 
for collisions with neutral hydrogen for simplicity, although the 
SCATE code includes tabulated coUisional cross-sections based 
on quantum-mechanical calculations for many species and tran- 
sitions (see AppendixfAli. No collisional broadening was applied 
to molecular line profiles. 

The effects of a background scattering continuum are 
demonstrated by analyzing the formation of fictitious Fe I lines 
at 3000 A in the atmosphere of the 3D model with [Fe/H] = 
-3.0. The left panel of Fig. |6] shows spatial averages of disk- 
center intensity profiles of a low-excitation (x - eV) transi- 
tion, where scattering opacity is treated as absorption (dashed 
line) and as coherent scattering (solid line). Line core brightness 
increases along with the continuum in the spatial average profile. 
The right panel of Fig. |6] shows the same profiles normalized 
to their individual average continuum intensity, as well as the 
line bisectors in the profile centers. Line depth increases when 
scattering is included, and the profile bisector shifts to slightly 
longer wavelengths. Note the characteristic C-shape of the bisec- 
tor, a consequence of spatial dominance of the upflowing gas in 
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Fig. 2. Spatial and temporal averages of the continuum flux distribution {F'^} as a function of wavelength, computed for time 
sequences of the 3D model atmospheres with metallicity -3.0 < [Fe/H] < 0.0 (upper left to lower right), treating scattering as 
absorption (dashed lines) and as coherent scattering (solid lines). The upper panel of each plot shows the deviation A {F^} of the 
coherent scattering cases from the continuum flux distribution where scattering is treated as absorption; vertical dotted lines indicate 
A - 3000 A, 4000 A and 5000 A where Fe line profiles were computed. Note the Balmer jump at 3647 A. 



the b right granules, which causes blueshifts in the emitted light 
(e.g.. lDravinslll982h . 

The depression of line core intensity in the normalized 
profile through continuum scattering stems from thermaliza- 
tion through line absorption. If line opacity is significant in 
continuum-forming layers, the thermalization depth of the ra- 
diation field moves outward into cooler parts of the atmosphere, 
since the joint photon destruction probability £y of continuum 
and line opacity is larger than the continuum photon destruction 
probabihty e^: 

ey = + -^^e' = 4^ ^ (7) 

line opacity is treated as absorption in our computations. The 
radiation temperature in the deeper parts of the line contribu- 
tion function decreases towards the lower local gas temperature, 
and the line core gains less brightness through scattering than 
the continuum. Normalization turns this disproportional growth 
into a deeper line profile. If line opacity is not significant in 
continuum-forming layers, line core brightness increases along 
with the continuum through the contribution from scattered pho- 
tons, and the intensity gain is mostly divided out when the line 
profile is normalized. 

Both magnitude and height-dependence of the line opacity 
thus determine the importance of continuum scattering for line 



formation: the ability of lines to thermalize the radiation field 
grows with x\'^ the opacity contribution of low-excitation lines 
of neutral atoms such as the Fel line in Fig.|6]is biased towards 
higher atmospheric layers (compare the distribution of line-to- 
continuum opacity rsAios Ix'^ the continuum optical surface 
T*^ = 1.0 and at t'^ =0.1 in Fig.|7]i, such species therefore have a 
weaker effect on the thermalization depth in continuum-forming 
layers. High-excitation lines exert larger influence, as the ex- 
ponential temperature-dependence of the Boltzmann excitation 
equilibrium moves their opacity contribution closer to t*^ = 1 .0. 

The growth of continuum-to-line core intensity ratios 
through scattering correlates with normalized continuum inten- 
sity across the steUar surface (left panel in Fig. [8]i: in the hot 
granules, lines gain less strength since line opacity is relatively 
small and continuum absorption relatively large in continuum- 
forming layers (left panel of Fig. Q and right panel of Fig. |5]l, 
and both continua and cores increase in brightness (left panels 
of Fig. |5] and Fig. [8]l. In the cool intergranular lanes, continua 
gain more intensity than the cores since line opacity plays a more 
important role in thermalizing the radiation field. 

The spatial dependence of line growth is translated into 
wavelength space as Doppler-shifts change sign between the up- 
flowing gas in granules and downflowing gas in intergranular 
lanes. The right panel of Fig. [8] shows the correlation between 
the growth of equivalent widths Wa through scattering and wave- 
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Fig. 3. Continuum photon destruction probabilities for arbitrary snapshots of the 3D model atmospheres with [Fe/H] = -3.0 (upper 
left panel), [Fe/H] = -2.0 (upper right panel), [Fe/H] = -1.0 (lower left panel), and [Fe/H] - 0.0 (lower right panel), averaged 
over surfaces of constant monochromatic continuum optical depth rjj for wavelengths X - 3000 A (solid lines), X - 4000 A (dashed 
lines) and A - 5000 A (dot-dashed lines). The monochromatic optical surfaces at rjj = 1 are marked by vertical dotted lines. 
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Fig. 4. Continuum surface intensities in the disk center at 3000 A computed for an arbitrary snapshot of the 3D model with [Fe/H] = 
-3.0, treating scattering as absorption (left panel) and as coherent scattering (right panel). Intensities are normalized to the average 
continuum intensity with coherent scattering. 



length shifts AAq of the individual line profiles in each column. 
We define AAq through the deviation of the intensity profile bi- 
sector at half-depth from the center wavelength Aq, 



Aq, 



(8) 



since the 3D velocity field may produce blended profiles with 
multiple minima in some columns. The stronger gains in line 
strength in the intergranular lanes are shifted towards larger 
wavelengths, causing the red wing of the normalized profile to 
appear relatively darker 

The flux profiles in Fig. |9] demonstrate that the eff'ects of 
scattering also become stronger towards the limb: continuum 
optical surfaces move outward into cooler layers with smaller 



continuum photon destruction probability (see Fig.O, which in- 
creases the brightness gain in the continuum and therefore pro- 
duces stronger lines (right panel of Fig.|9]l through the absorption 
mechanism that was discussed above. 



We quantify the impact of scattering on flux profiles by com- 
puting synthetic curves of growth, which are sampled with a set 
of typically 9 oscillator strengths and interpolated using cubic 
splines, treating scattering as absorption (dashed line in the left 
panel of Fig. [TOl i and as coherent scattering (solid line). The 
deeper profiles of the scattering case lead to larger equivalent 
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Fig. 5. Left: Correlation between the continuum intensity ratio /'(scat)//'(abs) and the normalized continuum intensities with scat- 
tering /'(scat)/ (/'(scat)) at 3000 A, computed for an arbitrary snapshot of the 3D model with [Fe/H] = -3.0. Right: Correlation 
between the continuum photon destruction probabilities e'^ at the local monochromatic continuum optical surface t'^ - \ and the 
normalized continuum intensities with scattering /'(scat)/ (/'(scat)) at 3000 A for the same model. 




-100 -50 50 100 -100 -50 50 100 

AX in [mA] AX in [mA] 

Fig. 6. Spatial averages of disk center intensity profiles (left panel) and normalized disk center intensity profiles with bisectors (right 
panel) of a fictitious low-excitation (x - ^ eV) Fe I line plotted as functions of wavelength shift A/l in mA with Aq - 3000 A, 
computed for an arbitrary snapshot of the 3D model with [Fe/H] - -3.0 including coherent continuum scattering (solid lines) and 
treating scattering as absorption (dashed lines). 




Fig. 7. Correlation between wavelength-integrated line-to-continuum opacity ratios x Ix'^ normalized continuum intensities 
/'(scat)/ (/'(scat)) at the continuum optical surface t' = 1 (left panel), and at continuum optical depth t' =0.1 (right panel) for the 
Fel line of Fig.|5](3D model with [Fe/H] = -3.0). 



widths Wx- The distance between the two curves at each W,\ de- 
fines a logarithmic 3D-3D abundance correctior0 

A log e = log e(3D, scat) - log e(3D, abs); (9) 

the solid line in the right panel of Fig. [TOl shows the result. As 
larger thermalization of the radiation field near line core wave- 

^ In spectroscopic notation, log 6(A) = log[Q(A'A/A'H) + 12.0, where 
A'a and A'h are the number densities of species A and hydrogen 



lengths produces increasingly deeper flux profiles, saturation in 
the curve of growth is delayed and stronger lines exhibit larger 
abundance corrections. Aloge also grows with excitation level 
X of the transition (dashed line in the right panel of Fig. [TOl l: 
larger line opacity in continuum-forming layers leads to abun- 
dance corrections of up to -0.36 dex at Wa - 60 mA. The ef- 
fect of continuum scattering on abundance measurements can 
easily be verified using an order-of-magnitude estimate, assum- 
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Fig. 8. Correlation between the growth of continuum-to-line core intensity ratios (Z*^// )scat/(^'^/^ )abs through scattering and nor- 
malized continuum intensities /'(scat)/ (/'(scat)) (left panel), as well as between the growth of equivalent widths W^(scat)/W^i(abs) 
through scattering and the line shift A/Iq (right panel) for the Fe I line of Fig.|6](3D model with [Fe/H] - -3.0). 
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Fig. 9. Spatial averages of flux profiles (left panel) and normalized flux profiles with bisectors (right panel) of the Fe I line of Fig. |6] 
(3D model with [Fe/H] = -3.0). 
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Fig. 10. Left: Curves of growth for the Fe I line of Fig.|9](3D model with [Fe/H] - -3.0) as a function of oscillator strength loggf, 
computed with scattering as absorption (dashed line) and as coherent scattering (solid line). Right: Corresponding 3D-3D scattering 
abundance corrections as a function of equivalent width Wa for the curves of growth shown in the left panel where = eV (solid 
line), and for a Fe I line with;if = 2 eV (dashed line). 



ing a rectangular spectral line with suflicient absorption opac- 
ity that the core flux is insensitive to continuum emission. 
Scattering then increases only the continuum flux F^, and the 
ratio of equivalent widths W,i(scat)/W^(abs) is independent of 
line width. For a linear curve of growth and an initial profile 
depth d = (I - F°/F'(abs)) = 0.5, a continuum flux gain 
F'(scat)/F'(abs) = 1.5 through scattering then translates into 
an abundance correction of the order of 0. 1 dex. 



7. 3D-3D scattering abundance corrections 

7.1. Curves of growth for Fe I and Fe II lines 

We compute scattering abundance corrections for fictitious Fe I 
and Fe II lines with excitation potentials x - and = 2 eV 
at 3000 A, 4000 A and 5000 A and for all 3D models to in- 
vestigate the dependence on ionization stage, excitation poten- 
tial, wavelength and metallicity. Figure [12] shows the results 
for the [Fe/H] = -3.0 model. As expected from the discus- 
sion in Sect. |6] we find the largest A log e for the strongest 
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high-excitation Fell lines at 3000 A, where thermalizing opac- 
ity becomes most important in continuum-forming layers near 
line core frequencies and desaturation has the strongest effect. 
Abundance corrections reach A log e =s -0.5 dex at Wa = 60 mA 
(dot-dashed line in the upper left panel of Fig.fTSI). The correc- 
tions become less severe for weaker lines, neutral species and 
lower excitation levels. At 4000 A (upper right panel), scatter- 
ing effects are still significant; strong high-excitation Fe II lines 
reach Aloge ~ -0.09 dex at Wa = 60 mA. At 5000 A (lower 
left panel), scattering is negligible with A log e ^ -0.02 dex at 
Wa = 60 mA, which is smaller than typical abundance measure- 
ment errors. 

The results for the [Fe/H] - -2.0 model and the [Fe/H] - 
-1.0 model (Fig.[T3]and Fig.[T4b exhibit similar behavior as scat- 
tering opacity is important for continuum formation (Fig. |2] and 
Fig. O. Abundance corrections are again largest for both metal- 
licities at 3000 A and for the strongest high-excitation Fe II lines, 
reaching Aloge x -0.5 dex ([Fe/H] = -2.0) and Aloge » 
-0.4dex ([Fe/H] = -1.0) at Wa = 60mA. The smaller effec- 
tive temperature of the [Fe/H] = -1.0 leads to relatively larger 
scattering effects compared to the other models (see Sect.|5]l. 

The situation changes at solar metallicity ([Fe/H] = 0.0, 
Fig. [TSl l: scattering is weak around the optical surface at all 
wavelengths (lower right panel of Fig. |3]l, leading to generally 
smaller abundance corrections. At 3000 A, A log e reaches only 
-0. 1 dex for the strongest high-excitation Fe II lines and only 
-0.03 dex at small equivalent widths. At 4000 A and 5000 A, 
scattering can be neglected. 

7.2. Curves of growth for molecular lines 

Spectral lines from molecules are important tools for abundance 
measurements of carbon, nitrogen and oxygen in metal-poor 
stars. Computing synthetic line profiles and curves of growth re- 
quires solving a set of equilibrium equations to obtain population 
numbers of the different atoms and molecules (see AppendixlAt. 
In late-type stellar atmospheres, simple molecules such as car- 
bon monoxide (CO), the hydrides CH, NH and OH and cyanide 
(CN) form in sufficiently cool layers. CH, OH and NH molecules 
have observable transitions in the blue and UV wavelength re- 
gions. We follow Collet et al. (2007) and compute synthetic 
curves of growth for fictitious molecular lines at 3150 A (OH), 
3360 A (NH) and 4360 A (CH). 

The resulting abundance corrections for the [Fe/H] - -3.0 
model are shown in the lower right panel of Fig. [12] We find the 
largest scattering effect for the OH line due to its short wave- 
length, reaching -0.13 dex for the strongest lines (dashed line 
in the lower right panel of Fig. [TSl l. NH lines at 3360 A ex- 
hibit a slightly smaller abundance correction (dotted line) with 
Alog e = -O.lOdex at Wa = 60mA; CH lines at 4360 A experi- 
ence only weak influence from continuum scattering (solid line). 
The magnitude of A log e for molecular lines is very similar to 
Fe I lines, as molecules form mostly in higher, cooler layers of 
the atmosphere. 

The abundance corrections of molecular lines at [Fe/H] - 
-2.0 are similarly large as for the most metal-poor model, but 
they become less severe at [Fe/H] = -1.0 where the OH cor- 
rection reduces to -0.09 dex at Wa = 60 A. Aloge is almost 
negligible at [Fe/H] = 0.0, the strongest OH features reach a 
correction of -0.03 dex. 

Note that we vary oscillator strengths instead of elemen- 
tal abundances to calculate the curve of growth. The resulting 
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Fig. 11. 3D-3D scattering abundance corrections as a function 
of equivalent width Wa for a fictitious OH line at 3150 A com- 
puted for 3D models with [Fe/H] = -3.0 using 4 opacity bins 
(solid line) and using 12 opacity bins (dashed line). 

correction A log e is therefore an approximation, as shifts in 
the molecular equilibria are not taken into account. Decreasing 
log e(0) by 0.13 dex in the [Fe/H] - -3.0 model and increasing 
log gf by the same amount leads to a deviation of the resulting 
equivalent widths of < 2 %, the systematic errors in the abun- 
dance corrections of Fig. [T2] through Fig.[T5]are therefore small. 

Computin g line pr o files for the OH molecule, 

[Gon zalez Her nandez et al.[ ([201 Oh recently reported on a 
dependence of their 3D oxygen abundances on the opacity 
treatment in their 3D model atmospheres of metal-poor dwarf 
stars. The shallower temperature gradient of their model with 
12 opacity bins significantly reduces the 3D- ID abundance 
corrections due to the strong temperature-dependence of 
molecular equilibrium populations. We therefore investigate the 
importance of the opacity treatment for the 3D-3D scattering 
effects of the OH 3150 line at [Fe/H] = -3.0. The resulting 
abundance corrections as a function of equivalent width for 
the 4 bin 3D model and the 12 bin 3D model (see Sect. [3]l are 
shown in Fig. [TT] The overall dependence of Aloge on the 
opacity treatment is very small, as the temperature-decreased 
OH population of the 12 bin model is compensated by increased 
abundance (or oscillator strength in our case) on the curve of 
growth to reach a given equivalent width. The differences in the 
atmospheric stratification become slightly more noticeable for 
stronger lines, which sample a wider height range, but the devi- 
ation between the predictions of A log e remains below 0.01 dex. 
The relative effects of scattering on abundance measurements 
should therefore be insensitive to the opacity treatment used for 
constructing the model atmosphere when comparing 3D line 
formation computations. 

8. 3D-1D abundance corrections 

Abundance analyses with ID hydrostatic model atmospheres are 
still commonplace in astrophysical research. We therefore com- 
pare the curves of growth derived from the 3D models with the 
results obtained for ID hydrostatic MARCS models. 

In order to establish the importance of the 3D structure and 
velocity fields, we first determine the scattering effects on a 
Fe I line with = eV that forms in a 3D atmosphere with 
zero gas velocities and in a ID hydrostatic MARCS model at- 
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Fig. 12. Upper left panel to lower left panel: 3D-3D scattering abundance corrections for fictitious Fel lines and Fell lines with 
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right panel: 3D-3D scattering abundance corrections for typical CH, NH and OH lines for the same model. 



0.00 
-0.10 
-0.20 
-0.30 
-0.40 
-0.50 
-0.60 

0.00 
-0.01 



Fe X=3000 A [Fe/H] = -2.0 „ Fe X=4000 A 



Fe I eV 
Fe I 2 eV 



Fe 
Fe 



eV 
2 eV 



30 40 
W, in [mA] 



Fe A=5000 A [Fe/H] = -2.0 



Fe I eV 
Fe I 2 eV 



Fe II eV 
Fe II 2 eV 



30 40 
W> in [mA] 



0.00 

u 

.E -0.02 
in 

o -0.04 

O 

ro 

% -0.06 
o 

i -0.08 

D 
(/) 

o" -0.10 

tr, 

-0.12 

o 

£ 0.00 



[Fe/H] = -2.0 



Fe I eV 
Fe I 2 eV 



Fe II eV 
Fe II 2 eV 



30 40 
Wi in [mA] 



CH NH OH 



[Fe/H] = -2.0 



-0.05 



-0.10 



CH 4360 A 
NH 3360 A 



30 40 
in [mA] 



Fig. 13. Same as Fig. [12] but computed for the 3D model with [Fe/H] - 



-2.0. 



12 



W. Hayek et al.: 3D LTE spectral line formation with scattering in red giant stars 



[dex; 


0.00 ! 


c 


-0.10 h 


tn 




n 
o 


-0.20 [ 


q' 

ili 




Ul 

CP 
O 


-0.30 r 


3C0t)- 


-0.40 r 




-0.50 [ 


oi 
o 


-0.60 ^ 



Fe A=3000 A [Fe/H] = -1.0 



Fe I eV 
Fe I 2 eV 



Fe II eV 
Fe II 2 eV 



30 40 
W> in [mA] 



Fe X=5000 A 



[Fe/H] = -1.0 „ 



Fe I eV 
Fe I 2 eV 



Fe 
Fe 



eV 
2 eV 



10 



20 



30 40 
W, in [mA] 



50 



60 



Fe A=4000 A 



S 0.00 

■D 

.E -0.02 
in 

D -0.04 
Q 
1^ 

'% -0.06 

O 

i -0.08 
o 

Q -0.10 
-0.12 



[Fe/H] = 



1.0 



Fe I eV 
Fe I 2 eV 



Fe II eV 
Fe II 2 eV 



30 40 
in [mA] 



CH NH OH 



[Fe/H] = -1.0 



-0.05 



CH 4360 A 
NH 3360 A 



20 



30 40 
Wi in [mA] 



50 



60 



Fig. 14. Same as Fig. [12] but computed for the 3D model with [Fe/H] - -1.0. 
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Fig. 15. Same as Fig. [12] but computed for the 3D model with [Fe/H] - 0.0. Note the different scaling of the vertical axes. 
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mosphere with [Fe/H] — -3.0. Missing broadening through 
gas motion is replaced with microturbulent broadening in these 
two cases, choosing ^ - l.Okms"', ^ = 2.0kms"' and ^ - 
3.0kms Oscillator strengths for the calculations with the ID 
model need to be increased by 0.7 dex to obtain approximately 
the same range of equivalent widths as in the 3D case. The ad- 
justment compensates for lower absorber population numbers, 
which again stem from the much shallower temperature gradient 
above the surface of the ID model compared to the 3D atmo- 
sphere; low-excit ation transitions at low metallicity are particu- 
larly affected (see Collet et al]|2007h . In the absence of Doppler- 
shifts, profiles are exactly symmetric (center column and right 
column of Fig.[T6]l. Resulting abundance corrections are shown 
in Fig.[T7] In the 3D static case, emission from granules and in- 
tergranular lanes is no longer separated in wavelength. At low 
microturbulence (^ = 1.0 km s"'), the line profiles saturate ear- 
lier compared to the 3D case with velocity field. Desaturation 
through scattering effects is reduced as the more weakly scatter- 
ing granules dominate the radiative flux, resulting in less depen- 
dence of A log e on line strength (solid line in the left panel of 
Fig.fTTl). Increasing microturbulent broadening delays saturation 
of line profiles in the scattering and absorption cases, extending 
the linear part of the curve of growth, and reduces thermaliza- 
tion as line opacity is distributed over a wider profile. Aloge 
is almost constant (dotted line and dashed line in the left panel 
of Fig. [TTI i. The weakest lines at Wa - 10 mA are indepen- 
dent from microturbulent broadening and yield approximately 
the same abundance corrections for all values of ^. 

The ID MARCS model exhibits similar scattering abundance 
corrections at the lowest equivalent widths with Aloge - 
-0.1 dex. They are more severe for stronger lines compared 
to the 3D static case, reaching -0.2 dex at Wa - 60 mA for 
^ - l.Okms"' (right panel of Fig. [TTll. Increasing microturbu- 
lence again desaturates both the absorption and scattering line 
profiles and reduces thermalization, resulting in a weaker de- 
pendence of A log e on line strength. 

We repeat the profile computation for the same fictitious 
lines as in Sect. |7] using the ID MARCS models and compare 
3D-1D abundance corrections with coherent scattering and 
treating scattering as absorption. Microturbulent broadening was 
fixed at ^ = 2.0 km s"' to simplify the discussion. Contrary to the 
differential comparison of 3D curves of growth in the previous 
section, the absolute 3D-1D abundance corrections in each case 
are more sensitive to the stratification of the model atmospheres. 
However, a test using our 12 bin model with [Fe/H] = -3.0 
for the OH line showed that the deviation between the 3D-1D 
abundance corrections with scattering and treating scattering as 
absorption was insensitive to the temperature profile of the 3D 
model. The following discussion thus focusses on the relative 
effect of scattering. 

The upper left panel of Fig.fTSlshows the 3D- ID abundance 
corrections for the Fe I and Fe II lines at 3000 A for the models 
with [Fe/H] - -3.0, computed with coherent scattering (black 
lines) and treating scattering as absorption (gray lines). The over- 
all behavior of the scattering-as-absorption calculations follows 
the results of Collet et al. (2007), who used model atmospheres 
with comparable stellar parameters in their analysis, and will 
therefore not be further discussed. The 3D- ID abundance cor- 
rections with scattering are very similar to the 3D-1D absorp- 
tion results for all weak lines as the 3D-3D and ID- ID scat- 
tering effects reach similar strength (» -0.1 dex in either case, 
see the right panel of Fig. [10] and Fig. [TTT i. As Wa grows, the 
calculations predict stronger scattering corrections in 3D, which 



leads to a growing downward deviation from the scattering-as- 
absorption 3D-1D results for all species. High-excitation Fell 
lines exhibit the largest difference, requiring an adjustment of the 
3D- ID scattering-as-absorption abundance corrections by up to 
0.35 dex. It is important to keep in mind that the effect of scat- 
tering on abundance corrections in ID depends on the choice of 
microturbulent broadening (right panel of Fig. [TTl i: choosing a 
smaller ^ thus reduces the deviation between the 3D- ID abun- 
dance corrections. However, the general influence of microtur- 
bulence on 3D-1D corrections through the ID curve of growth 
is often larger than that (> 0.1 dex for low-excitation Fel lines 
with Wa - 60 mA when f is in creased by 0.5 km s"'; see Fig. 5 
and Fig. 6 in Collet et al.ll2007h . 

At longer wavelengths, the importance of scattering de- 
creases, and the 3D- ID abundance corrections with scattering 
come into closer agreement with the scattering-as-absorption 
calculations for all species and line strengths (upper right panel 
and lower left panel of Fig. [TsT i: the deviations reach 0.05 dex 
at 4000 A for the strongest Fell lines and almost vanish at 
5000 A. The same is observed for CH, NH and OH molecular 
lines, with OH exhibiting the largest deviation of 0.05 dex be- 
tween the 3D- ID corrections when scattering is included for the 
strongest lines, while the differences are negligible in the case of 
CH lines. The treatment of molecular lines in a 3D- ID com- 
parison is more complex due to the strong non-linearities that 
stem from the temperature sensitivity of the equilibrium popula- 
tions and their dependence on the chemical composition of the 
model, which may change dramatically at the lowest metallic- 
ities (Collet et al. 2007). We fix the abundances of carbon, ni- 
trogen and oxygen to the scaled solar composition; this leads to 
an offset of the absolute 3D - ID coiTections with respect to the 
results of ICollet et"d] (l2007h . but only weakly affects the differ- 
ential comparison of scattering effects. The 3D-1D abundance 
corrections for molecules should therefore not be directly ap- 
plied to ID abundance analyses. 

At metaflicities [Fe/H] = -2.0 and [Fe/H] = -1.0, the pre- 
dicted 3D- ID abundance corrections generally decrease due to 
the growing similarity between the temperature stratifications of 
the 3D models an d the ID models (for further discussion see 
ICoUet et alJl2007l) . but the deviations of the scattering cases are 
very similar to the calculations with the [Fe/H] = -3.0 mod- 
els (Fig. [19] and Fig. l20t . as scattering is still important at the 
shortest wavelengths and for the strongest lines (see Fig. [T3] and 
Fig.[T4]i. At solar metallicity, the 3D-1D abundance corrections 
for the strongest high-excitation Fe II lines at 3000 A deviate by 
less than 0. 1 dex when scattering is included (upper left panel of 
Fig.l2ni. and they are negligible at longer wavelengths. 

The effects of the 3D temperature stratification and 3D struc- 
ture dominate the 3D-1D abundance corrections at all metaflic- 
ities for all Fe lines with A > 4000 A and for molecular lines, 
but scattering radiative transfer nevertheless leads to significant 
deviations for the strongest Fe lines at 3000 A. It is not clear 
whether this removes or adds trends with wavelength in the 
3D and ID line-by-line abundance measurements for a given 
species, since the magnitude of the deviation depends on line 
strength and on the microturbulence parameter which is needed 
to broaden ID line profiles; detailed stellar abundance measure- 
ments will be needed for further insight. 

9. Conclusions 

We investigate the effects of a scattering background continuum 
on spectral line formation based on 3D hydrodynamical and ID 
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Fig. 16. Spatial averages of flux profiles as functions of Doppler shift Av from the line center for the Fel line of Fig.|9l computed 
treating scattering as absorption (dashed lines) and as coherent scattering (solid lines) using the 3D model with [Fe/H] = -3.0 
(left column), the 3D model with [Fe/H] - -3.0 and with velocity fields set to zero (3D static model, center column) and the ID 
hydrostatic MARCS model with [Fe/H] = -3.0 (right column). The bottom row shows normalized profiles. The 3D static model and 
the ID model assume microturbulent broadening with ^ = 2.0 kms '; line strengths of the ID calculations needed to be increased 
by +0.7 dex to obtain similar equivalent widths (see text). 




Fig. 17. Abundance corrections for the same Fel line as shown in Fig.|9] but computed for the 3D model with [Fe/H] - -3.0 and 
with velocity fields set to zero (left panel) and for the corresponding ID hydrostatic MARCS model with [Fe/H] = -3.0 (right panel), 
assuming microturbulent broadening with ^ - l.Okms ' (solid lines), ^ - 2.0kms ' (dotted lines), and ^ = 3.0 kms ' (dashed 
lines). 



hydrostatic model atmospheres of red giant stars with different 
metallicity. Continuum scattering is treated in the coherent and 
isotropic approximation; opacities assume LTE. We compute ra- 
diative transfer using the SCATE line formation code to determine 
synthetic continuum flux levels and flux profiles of fictitious Fe I 
and Fe II lines between 3000 A < A < 5000 A, as well as selected 
fictitious CH, NH and OH molecular transitions. 

Rayleigh scattering contributes significant opacity in the 
blue and UV wavelength bands, allowing photons to escape from 
larger atmospheric depths compared to the case of treating scat- 
tering as absorption. As a consequence, the continuum surface 
flux increases above 2000 A and below 5000 A, while scatter- 
ing is negligible in the infrared and at even shorter wavelengths, 
where photoionization processes thermalize the radiation field. 
In the 3D models, the strongest relative brightness gain is ob- 
served in intergranular lanes, where scattering opacity is more 



important and photon destruction probabilities near the optical 
surface are smaller compared with the hot granules. 

Continuum scattering affects the strength of LTE lines, de- 
pending on ionization stage, excitation potential, wavelength 
and metallicity. If a line contributes significant absorption opac- 
ity in continuum-forming layers, the thermalization depth near 
line core frequencies moves outward into cooler layers and the 
core gains less brightness through scattering than the contin- 
uum. Normalization of the flux profiles translates this dispro- 
portional intensity gain into a deepening of the line. We find 
the largest effects at the lowest metallicity, the shortest wave- 
lengths, for singly ionized lines and for higher excitation levels, 
where continuum scattering is strongest and line opacity is sig- 
nificant in continuum-forming layers. The temperature contrast 
across the surface of 3D model atmospheres results in diff'eren- 
tial line growth, as intergranular lanes experience a stronger scat- 
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tering effect than granules. Doppler-shifts reverse sign between 
granules and intergranular lanes, translating spatial variation into 
wavelength space. The red wing of the lines appears thus more 
strongly deepened than the blue wing, and the line center shifts 
slightly towards longer wavelengths. 

We quantify the importance of continuum scattering in 3D 
by comparing synthetic curves of growth using scattering ra- 
diative transfer and treating scattering as absorption. Increasing 
line strength through scattering desaturates line profiles, lead- 
ing to a growing deviation between the curves of growth for 
stronger lines. Transformed into a 3D-3D abundance correction 
for given equivalent width, we find up to x -0.5 dex deviation 
forFellwith;^' = 2 eV at 3000 A and [Fe/H] = -3.0. At 4000 A, 
the effects of scattering on line abundances are much smaller 
and only significant for the strongest lines; at 5000 A, scattering 
can be neglected. Transitions of the CH, NH and OH molecules 
behave in an analogous way to neutral low-excitation atomic 
lines; the strongest 3D-3D abundance corrections appear at the 
shortest wavelengths. Scattering effects weaken towards higher 
metallicity as continuous absorption opacity becomes increas- 
ingly important; 3D-3D abundance corrections at solar metal- 
licity ([Fe/H] = 0.0) only reach -0.1 dex for the strongest high- 
excitation Fe II lines at 3000 A, while they are practically negli- 
gible for weak lines and at longer wavelengths. 

The importance of velocity fields for line formation in the 
given wavelength range is tested by comparing the results with a 
3D atmosphere with zero velocities and artificial microturbulent 
line broadening. In the absence of Doppler-shifts through the 
granulation flow, the line profiles become symmetric and satu- 
rate earlier for smaller micro turbulence. The desaturating effect 
of scattering is weaker as granules dominate the emission, where 
photon destruction probabilities are larger, resulting in smaller 
abundance corrections. We repeat the curve of growth computa- 
tion for ID hydrostatic MARCS models with the same stellar pa- 
rameters as the 3D models and find that scattering leads to a sig- 
nificant downward adjustment of the 3D-1D abundance correc- 
tions for the strongest Fe lines at the shortest wavelengths com- 
pared to 3D-1D corrections that were computed treating scat- 
tering as absorption; the deviation reaches -0.35 dex for high- 
excitation Fell lines at 3000 A, [Fe/H] = -3.0 and microtur- 
bulence ^ = 2.0km s"'. Weak lines are much less affected, as 
3D-3D and ID- ID scattering abundance corrections are very 
similar. The scattering corrections for molecules are small com- 
pared to the 3D effects as they remain below 0.1 dex in all cases; 
our absolute 3D- 1D corrections are offset from the results of 
ICollet etal.l(l2007h as we use fixed molecular equilibrium popu- 
lations. 

Chemical abundance analyses of metal-poor giant stars that 
include spectral lines in the UV and blue regions should take 
background continuum scattering into account. Convective ve- 
locity fields change the profile shapes, which shift towards 
slightly longer wavelengths and produce larger scattering abun- 
dance corrections for stronger lines compared to ID hydrostatic 
models; it is therefore important to conduct abundance analyses 
based on 3D hydrodynamical model atmospheres. 

Our investigation assumes coherent isotropic scattering, ne- 
glecting redistribution of radiation through thermal Doppler- 
shifts, which should only have a small effect on the synthetic 
line profiles. The LTE approximation is a more severe limita- 
tion: non-LTE effects, such as photoionization, are known to be 
important in metal-poor stars and need to be considered for a 
quantitative an alysis of stellar hne profiles (for a discussion, see 
lAsnlundll2005h . 
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Appendix A: The SCATE code 

The SCATE spectral line formation code for stellar atmospheres computes 3D ra- 
diative transfer for blended spectral lines and a background continuum with co- 
herent isotropic scattering. Local thermodynamic equilibrium (LTE) is assumed 
to compute level populations for lines and continuous opacity sources, using 
Saha-Boltzmann equilibrium calculations. A lookup-table is used for monochi'o- 
matic continuum opacities, continuum photon destruction probabilities, Saha 
ionization equilibria, molecular equilibria (for hydrogen, carbon, nitrogen and 
oxygen) and partition functions to speed up the 3D line formation computations, 
in particular when a time-series of profiles is required. 



AA. Numerical solution of the radiative transfer equation 

Computation of a line profile is divided into two steps: a short-characteristics- 
based radiati ve transfe r solver with Gauss -Seidel-type convergence acceleration 
(Trujillo Bueno & Fabiani B endicholll995 ;^ computes the monochromatic radi- 
ation field at each frequency point v across the profile to obtain the contin- 
uum source function with coherent scattering. The Gauss-Seidel method uses an 
upper/lower triangular approximate A* operator which computes source func- 
tion corrections during the formal solution and yields fast convergence; see 
[Havek et al. (2010) for a detailed description of the implementation. The solver 
takes the anisotropy of the combined continuum and line source function in the 
observer's frame into account, which is due to the angle-frequency coupling of 
line opacity through Doppler-shifts. Linear or higher-order local inteipolation is 
available for translating opacities, source functions and upwind intensities onto 
the characteristics grid. The discretized source function integral is computed us- 
ing linear inteipolation or second-order inteipolation, which is required at high 
optical depths to correctly reproduce the dift'usion approximation. 

A second solver delivers angle-resolved surface intensities an d surface 
fluxes for obtaining the line profiles. SCATE offers the differential iFeautrieJ 
Q963) scheme as well as an integral sch eme which is based on the formal solu- 
tion of the transfer equation I Olson & Kunasdl 19871) . The integral scheme com- 
putes the source function integral using the same interpolation methods as the 
scattering solver The radiative transfer equation is solved on long characteris- 
tics that begin in the top layer of the hydrodynamical mesh and span across the 
simulation domain, until they reach the optically thick dift'usion region. Before 
each calculation of surface intensities, the 3D model atmosphere is tilted into the 
ray direction using local cubic interpolation. 

Both solution steps are parallelized over angle using OpenMP directives for 
shai'ed-memory ai'chitectures, allowing for fine resolution of the angle-frequency 
couphng in the scattering computation and for integrating surface intensities over 
a large number of rays when radiative fluxes are needed. A modern 8-core cluster 
node delivers a line profile with scattering (spatial resolution 120 X 120 X 230 
points, frequency resolution 40 points, angular resolution for the scattering 
solver 24 ray directions, angular resolution for flux profile computation 16 ray 
directions) in » 15min per atmosphere snapshot, while the radiative transfer 
problem without scattering is solved within < 1 min. 

SCATE uses partial grid refinement on the vertical axis to improve the resolu- 
tion of continuum optical surfaces. Additional layers are automatically inserted 
into the 3D model where vertical optical depth steps are large using cubic spline 
interpolation. The vertical axis is truncated when the optically thick dift'usion 
region is reached and the radiation field is entirely thermal. 
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Fig. 18. Upper left panel to lower left panel: 3D- ID abundance corrections for fictitious Fel lines and Fell lines with excitation 
potential = OeV and;\^ = 2eV at 3000 A, 4000 A and 5000 A, computed for the 3D and ID models with [Fe/H] = -3.0, including 
coherent scattering (black) and treating scattering as absorption (gray). The ID calculations assume microturbulent broadening with 
^ = 2.0 km s' . Lower right panel: Abundance corrections for typical CH, NH and OH lines for the same models. Note that feedback 
of the 3D abundance corrections on the molecular equilibrium is not taken into account (see text). 
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Fig. 20. Same as Fig. [18] but computed for the models with [Fe/H] - -1.0. 
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Fig.21. Same as Fig. [18] but computed for the models with [Fe/H] - 0.0. 
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A.2. Line opacities in tine LTE approximation 

The hne opacity at frequency v' in the local rest frame of a gas parcel in the 
stellar atmosphere is computed using the expression 



Xv' 



ne 

/lu«l 1 - 



e tT \ ifilv' 



(A.l) 



with the electron charge e, the electron mass m^, the speed of light c, the Planck 
constant h, the Boltzmann constant k and the temperature T. Literature values 
of the oscillator strength /j„ are usually combined with the statistical weight g\ 
of the lower level of the transition, which enters the opacity through Boltzmann 
excitation factors for the level population density «i . The code computes Saha- 
Boltzmann equilibria with tabulated partition functions that were calculated us- 
ing the NIST database (Ralchenko et alJ20 09'). Molecule formation is taken into 
account for 12 important species by computing LTE equilibrium populations 
(see Table lA. n . The term in brackets in Eq. )A. It corrects the line opacity for 
stimulated emission, assuming atomic level populations in LTE. 

The line profile if/ is given by a Voigt function around the laboratory wave- 
length yj) of the transition. Doppler-shifts through macroscopic velocity fields u 
in the observer's frame are included in the non-relativistic approximation, which 
yields the line profile 



i//(v-vo) = il/\v-v[ 



(A.2) 



in the observer's frame, for a ray in direction n and gas velocity u. The vectoriz- 
able recipe of Hui et al. tl978) is used to evaluate the Voigt profile numerically. 

Line profiles include Doppler broadening to account for shifts through ran- 
dom thermal gas motion. Microturbulent broadening may be included in the 
computations, but it is generally unnecessary for line formation with 3D hy- 
drodynamical atmospheres (see the discussion in Asplund et al. 2000b). 

Natural line broadening through radiative de-excitation is set by the profile 
width parameter E""' using laboratory data if available, or it may be estimated 
using the expression 



fiu. 



(A.3) 



which assumes a two level atom; A is the transition wavelength and is the 
statistical weight of the upper level. 

Pressure broadening through collisions of neutral atoms and ions with 
neutral hydrogen atoms is taken into account in the impact approxima- 
tion. Profile widths are based on quantum mechanical calculations by 
lAnstee & O'Mara (1995); Barklem & O'Mai'a (1997); Barklem et al. 1I9M 
boOOa) ; iBarklem & O'Mara >2001i) ; iBarklem & Aspelund-Johansson 12001 . 
CoUisional cross-sections are estimated using the power law 



o-(v) = o-(vo) — 



(A.4) 



where v is the relative velocity of the particle and the perturbing HI atom, (t(i'o) 
and a are tabulated, and I'o is a reference velocity. The width parameter of 
the Lorentzian broadening profile is obtained by integrating over a Maxwell- 
Boltzmann distribution for i', which yields 



•yHi = 2|-| E|— ^|vo-(vo)|-| nm, (A.5) 

with the gamma function F, the number density «hi of neutral hydrogen, and 
the average velocity v = ^J&kT/nfi, where /.( is the reduced mass of the parti- 
cle and the perturbing H I atom. The contribution from He I atoms to pressure 
broadening is estimated through 



/™«r"'|10 + CHeI 



»Hel \ 



(A.6) 



using an approximate scaling factor chcI ~ 0.41. In the absence of quantum 
mechanical data, the classical lUnsold 1 19551) approximation is used in the for- 
mulation found in Gray (2003): 



log r" ' = 20.0 + 0.4 log Ce + log P - 0.7 log T, 



(A.7) 



for quantities given in cgs units; P is the gas pressure. The interaction constant 
C(, is approximated through the expression 



Ce = 0.3 ■ 10" 



[xr 



(A.8) 



with the ionization potentials Xu" ™d x'"" from the upper and lower level of the 
transition, both given in eV. Eor the more c omplicated case of hydrogen lines, 
broadening recipes of lBarklem et alj j2000bl) are available. 



Table A.l. Molecular data included in the equilibrium cal- 
culations; partition functions were computed using the NIST 
database (Ralchenko et al. 2009). 



Molecule 


Dissociation energy 


H2 


Sauval & Tatum (1984) 




Sauval&Tatum(1984) 


H2O 


Ruscic et al. (2002) 


CH 


^uval& Tatum (1984) 


C2 


Urdahl etal. (1991) 


CN 


Huaneet al. (1992) 


CO 


mdelsberg etal. (1987) 


OH 


Ruscic et al. (2002) 


O2 


_Sauval& Tatum (1984) 


NH 


Marauette et al. (1988) 


N2 


~auval& Tatum (1984) 


NO 


Sauval & Tatum (1984) 



